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Abstract 

Energy dependence of heavy quarkonia production in hadron-nucleus collisions is 
studied in the framework of the Glauber-Gribov theory. We emphasize a change in 
the space-time picture of heavy-quark state production on nuclei with energy. Lon- 
gitudinally ordered scattering of a heavy-quark system takes place at low energies, 
while with increasing energy it transforms to a coherent scattering of projectile par- 
tons on the nuclear target. The characteristic energy scale for this transition depends 
on masses and rapidities of produced particles. For J/ip, produced in the central 
rapidity region, the transition happens at RHIC energies. The parameter-free cal- 
culation of J /ip in dAu collisions is in good agreement with recent RHIC data. We 
use distributions of gluons in nuclei to predict suppression of heavy quarkonia at 
LHC. 

Key words: J/tp absorption in nuclear matter, nuclear effects in hadron-nucleus 
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1 Introduction 



The heavy-ion programme at RHIC (BNL) and LHC (CERN) aims at discov- 
ering features of a possible new state of deconfined QCD matter anticipated 
to form in nucleus-nucleus collisions. An important signal of the quark-gluon 
plasma (QGP) formation would be a suppression of charmonium yield [1] pro- 
duced in these collisions. A proper baseline for the discovery of this effect is 
charmonium production in hadron-nucleus collisions, where the QGP is absent 
and only cold nuclear matter effects are present [2|3] . 
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Nuclear effects in hadron-nucleus collisions are usually discussed in terms of 
the power-law parameterization 



d°feA _ dg-fejy A a(x F ) /-.N 



where a% A (&h N ) is the inclusive cross section of particle a off a nucleus (nu- 
cleon). The function o,(xf) characterizes nuclear effects at different longitu- 
dinal momentum fractions of the produced particle, xf- For J/ip production, 
measurements show a decrease of a from 0.95 at xp ~ to values ~ 0.75 
at %f — 0.8 [4"|5f6] thus indicating an increase of absorption as xp increases. 
Also, data over a large range of energies show an approximate scaling of a 
with xp [415116] rather than a scaling with X2 (fraction of the total momentum 
carried by a parton from a nucleus) expected from QCD factorization. 

Recently, the substantial decrease of nuclear absorption in J/ip production 
in deuteron-gold (dAu) collisions at RHIC energy, \/s = 200 GeV, where 
Cabs ~ 1 — 2 mb [7] , compared to a abs ~ 4 mb j8] measured in proton- lead (pPb) 
collisions at SPS energy, y/s = 17.3 GeV has attracted a lot of attention. This 
corresponds to a value of a consistent with 1 at Xp = 0. It was widely believed 
that absorptive effects would increase or, at least, remain constant with rising 
collision energy [9|10llllj . In the model of [11] (Sec. 4.1), e.g., the postulated 
growth of a a b s with energy is motivated by a growth of charmonium-nucleon 
center-of-mass energy, reflecting the rapid growth of partons carrying small 
momentum fractions in the nucleon. 

An equally important implication resides in the fact that the RHIC data does 
not scale in xp: whereas the lower-energy data points display a flat behavior 
at small xp, the new points delineate a steep tilt. This novel feature seems 
also to be hard to reproduce in models describing the energy dependence of a 

mm- 

Such a behavior of a(xp) allows for a natural explanation in the Glauber- 
Gribov theory of multi-particle production on nuclei [H]. At very high energies 
Abramovsky-Gribov-Kancheli (AGK) cutting rules [15] lead to a cancellation 
of the Glauber-type diagrams in the central rapidity region, i.e. for Xp ~ 0, and 
only so-called "enhanced" diagrams [TBlfTT] . corresponding to multi-Pomeron 
interactions, contribute to a difference of a from unity. For light quarks, this 
coherent hadroproduction sets in at a typical energy scale E ~ m N fiRA, 
where Ra is the radius of the nucleus, is the mass of a nucleon and fi is 
a typical hadronic scale of the order of ~ 1 GeV. For heavy quark states, the 
mass M Q q of the heavy system introduces a new scale 
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where x + = \{Jx 2 F + AMq Q / s + xf) is the longitudinal momentum fraction 
of the heavy system. 

The AGK cutting rules are violated at larger values of xp and at low energies; 
the first effect is interpreted as conservation of energy- momentum [18]. The 
latter effect is in turn related to a change of the space-time picture of the 
interaction [19]. At energies below sm longitudinally ordered rescatterings of 
the heavy system take place. In this situation we can unambiguously define the 
production point of the heavy system and, in turn, the distance it has to travel 
through the surrounding nuclear matter. This leads naturally to the notion of 
an absorptive cross section. At s > sm the heavy state in the projectile, which 
also includes light degrees of freedom, scatters coherently off the nucleons of a 
nucleus, and the conventional treatment of nuclear absorption is not adequate. 

In the central rapidity region, the values of sm for J / ip are within the RHIC 
energy range. Accordingly, the effects of shadowing of nuclear partons become 
important and can be calculated using the Glauber-Gribov theory of nuclear 
structure functions in the region of x 2 < (mjyi?^) -1 . 



2 Model description and comparison to data 



Consider first production of heavy onia in the "low" energy regime, s < sm- 
It was shown in [TH] that the contribution of all diagrams with intermediate 
heavy-quark state to the total cross section is canceled in this energy region. 
However, the different s-channel discontinuities (cuttings) of these diagrams 
are different from zero. Consider, as an example, the cuttings shown in Fig. [TJ 
Their contributions to the total cross section are equal in magnitude and have 
opposite signs. Their contributions to the inclusive cross section would also 
cancel each other if the QQ final states in Fig. [T](a) and Fig. [H(b) would have 
the same distribution in xf- However, the elastic rescattering in Fig. [T^a) does 
not change the momentum, while in the case of inelastic interaction of the QQ 
system, Fig. d](b), it can loose some of its momentum and (or) transform into 
another state which is weakly coupled to the observed particle a (absorption). 
We parameterize the probability to produce these states in a single rescattering 
by a parameter (1 — e). 

Thus at energies s < sm the absorption is determined by 
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a) 



b) 



Fig. 1. (a) Elastic and (b) inelastic scattering of a QQ system. The former does not 
have to be the leading one. 



The function f% N (x+) = a^ N F^(x + ), where Ff denotes the unmodified distri- 
bution of produced particles, and £(x+) = (1 — e) + ex+ determines the x + 
dependence of the strength of the shadowing. Finally, 7^(6) = /f^ dz pA.(b, z) 
is the nuclear thickness function normalized to A. With £,(xp = 0) = (1 — e) 
one recovers the well-known Glauber formula [20|2T] with a a i, s = (1 — e) <j q q- 
As shown in Ref. [19], Eq. (j3J), with the inclusion of shadowing effects (see 
below), gives a good description of experimental data on charmonium produc- 
tion in pA collisions at E L ab ^ 800 GeV/c with a Q Q = 20 mb and e = 0.75. 
This corresponds to an absorption cross section of a a bs = 5 mb. Note, that 
a Q Q is rather large, indicating that the cc pair is produced in the color octet 
state rather than in the colorless state. It can also be viewed as a DD (D*D*) 
system. 

Equation (jHJ) is not applicable at asymptotic energies as the assumption of 
longitudinal ordering is only valid at s < sm- For energies higher than sm 
the expression will change due to the correct treatment of coherence effects 
according to [9] 



which is similar to the energy-momentum conservation effect for light quarks 
[T9] . In the model proposed in [9], a Q Q is equal to the total cross section of 
the QQ — N process, and is not proportional to x+. 

We would like to point out that this leads to an unnatural behavior at high 
energies due to the smallness of the Pomeron vertex and to an effective double- 
counting of nuclear effects, and propose an alternative procedure. If one con- 
siders non-enhanced Glauber-type diagrams, then the effective cross section 
varies as a Q Q ~ x+, thus satisfying the AGK cancellation. The suppression is 
concentrated at much higher xf for QQ production than for the light hadrons 
because of the large mass of the QQ system. It was shown in Ref. [19J that 
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Fig. 2. J/ip suppression as a function of centrality at different rapidities at RHIC. 
Data are taken from |7|31j . 

at xp ~ 1 the second rescatterings in the low and high-energy limits should 
coincide. This means that a Q Q ~ ex\a Q Q. Experiment on J/ip production in 
dAu collisions at RHIC [7] was performed in the central rapidity region, where 
x + varies from 0.025 to 0.05 and a Q g is, therefore, very small. This suggests 
an analysis of J/ip suppression in dAu collisions at RHIC energies taking into 
account enhanced diagrams only. A similar approach, albeit with a simpler 
parameterization of nuclear shadowing, has been considered in [22] . 

We have studied gluon shadowing in Ref. [23|241l25 j, where a model for 7* A 
collisions was considered within the Glauber-Gribov theory [H] including 



enhanced diagrams or, in other words, interactions among Pomerons. Sum- 
ming up an arbitrary number of Pomeron tree diagrams as in the generalized 
Schwimmer model [26] one obtains the following expression for the total cross 
section of a 7* A collision 



a s ^ h A (x } Q 2 ,b 
where 



Act, 



7* N 
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Fig. 3. Rapidity dependence of J/tp suppression for minimum bias dAu collisions at 
RHIC and predictions for pPb collisions at LHC. Data are taken from |7|31j . 

with xjp = 0.1, where shadowing is expected to disappear. Here F 2 (x, Q 2 ) is the 
structure function for a nucleon, F^(xp, Q 2 , (3) is the t-integrated diffractive 
structure function of the nucleon, B{xjp) is the t-slope of the diffractive distri- 

Xp. Equation ((Sj) 
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N 



bution, and Fji(t') is the nuclear form factor where t' 
determines shadowing for quarks (anti-quarks) in nuclei. For gluons the same 
expressions were used with substitutions F 2 £> (xp,Q 2 , 0) — > Fjp (xjp , Q 2 , (3) , 
F 2 (x,Q 2 ) — > xg (x,Q 2 ), indicating gluon distributions in the Pomeron, mea- 
sured in diffractive deep inelastic scattering (DDIS), and in the proton, re- 
spectively. Gluon distribution of the nucleon is taken from CTEQ6M param- 
eterization [27] . We take information on the diffractive gluon distribution and 
Pomeron parameters from recent HERA measurements [28], where two inde- 
pendent fits of the gluon diffractive distribution function, called FIT A and 
FIT B, represent the overall uncertainty of extracting this information from 
the measurements. We will show results only from the latter fit, denoted as 
GGB, as it is closer to, as yet, preliminary combined fits where di-jet produc- 
tion has been included [29J. This model has previously been used to calculate 
quark shadowing in nucleus-nucleus interactions in Ref. 



PHENIX collaboration has measured the nuclear modification factor (NMF) 
of J/ij) production in dAu collisions at RHIC as a function of centrality and 
rapidity in [7J and most recently in [3T| (with a more up-to-date pp reference). 
We define the centrality dependent NMF as 



RdAu ({Neon)) 



Njj u ((N coll )) 
(Ncou) x N'. 



inv 
VP 



(6) 
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where the average number of nucleon-nucleon collisions (N co u) is obtained from 
the Glauber model for a given centrality. The results of calculations based on 
Eq. (j^D are shown in Fig. [2] for the NMF, given by Eq. at backward, mid- 
and forward rapidity. Since the model of gluon shadowing does not include 
anti-shadowing effects, the result is quite trivial in the backward hemisphere, 
although not inconsistent with the data. Anti-shadowing is assumed to be a 
10% effect. At rapidity y = and y = 1.8 the consistency with experimental 
data is quite good. The rapidity dependence of nuclear modification factor 
RdAu for minimum bias dAu collisions at RHIC [7f3"l~] and predictions for pPb 
collisions at LHC, y/s = 5.5 TeV, are presented in Fig. [3J At mid-rapidity, 
gluon shadowing at LHC is a 40% effect, being barely significant, ~ 10%, 
at RHIC. This fact is important for the calculation of charmonium yield in 
nucleus-nucleus collisions at these energies. 

The current PHENIX data [31] indicate that at forward rapidity shadowing 
alone cannot be accounted for the full drop of the NMF of J/ip. Therefore, the 
dash-dotted curves in Fig.[2]and Fig.[3]depict calculations including both gluon 
shadowing and a model for energy-momentum conservation to be presented in 
the next section. The latter effect is shown to be relevant already at y = 1.8 
at RHIC. 



3 Energy dependence of a(ip) 

Based on the previous discussion we will now formulate a model for J/ip 
production in hadron-nucleus collisions at all energies. Figure H] shows exper- 
imental points for a in Eq. |T| as a function of both xp and x<i [4,5.7]. In 
order to calculate J/ip suppression at different energies and for all values of 
xf it is necessary to make a modification of Eqs. (EJHJ) to describe a transition 
from low-energy to high-energy regime, which happens at x c 2 ~ im^RA) 1 ■ 
The term ex\ in the expression for (or (t Q q(x + )) provides a smooth 

transition between the two regions. On the other hand, the term (t Q q(1 — e) 
describing the absorption should be modified in such a way that it should tend 
to zero in the high-energy region, being substituted in this region by gluon 
shadowing. To fulfill these requirements we introduce an extra multiplier 



for this term. This procedure corresponds to a suppression of Glauber-type, 
non-enhanced diagrams at high-energies, and has the correct high-energy be- 
havior, i.e., it satisfies the AGK cutting rules. 

The dashed curve in Fig. H] at y/s = 39 GeV has been calculated this way with 




(7) 
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Fig. 4. a dependence in pA collisions vs. x 2 (left) and xp (right) for y/s = 39 — 5500 
GeV. Data are taken from [4|5|7j . 

e = 0.93 and o q q = 45 mb. Following the prescription of Eq. for the high- 
energy regime, solid and dotted curves are calculations of gluon shadowing 
and energy conservation effect for RHIC and LHC energies, respectively. The 
calculations for RHIC have been performed for all < xp < 0.8, while for 
LHC we have, for illustrative purposes, only made calculations for the central 
rapidity region, \y\ < 4. Although the coverage in xf is small at this energy, 
the structure function of the nucleus is probed down to X2 ~ 10~ 5 . 

The model perfectly reproduces Xp, x 2 and energy dependence of all the ex- 
perimental data. Going from low to high energies we observe a breaking of 
xp scaling in the central rapidity region. The origin of this is that absorption 
effects, related to longitudinally ordered rescatterings, die out while shadow- 
ing slowly appears. This is most clearly seen in the a vs. x 2 plot. The form 
of the curve for a(xp) is reinstated at Xp > 0.25, although shadowing leads 
to a stronger overall suppression. Additionally, scaling in x 2 is predicted to 
appear for J/ip in the common kinematical window of RHIC and LHC, i.e. 
for 10~ 3 < x 2 < 0.05. This is a novel feature in heavy-ion experiments and 
would imply the validity of the factorization theorem in hadronic processes at 
ultra-relativistic energies. 

Concluding, we have argued that recent data on J/ip production at RHIC 
imply a profound change of the space-time picture of charmonium produc- 
tion in hadron-nucleus collisions and described the experimental data at mid-, 
forward and backward rapidities in terms of nuclear shadowing. Already at 
y > 1.7 we obtain a quite strong inflence of energy-momentum conservation, 
in accordance with most recent data from experiment. Furthermore, we pre- 
sented a model to describe the energy dependence of these features at xp 
larger or equal zero. The agreement with available data is very satisfactory. 
Nuclear effects at < are out of the scope of this paper, and require a 
separate study [32] • These findings confirm the appearance of shadowing ef- 
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fects in light particle production in dAu collisions, and will also have a great 
impact on models for nucleus- nucleus collisions both at RHIC and LHC. 
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